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It is shown that experimentally obtained isotherms of ad- 
sorption on solid substrates may be completely reconciled 
with Lifshitz theory when thermal fluctuations are taken into 
account. This is achieved within the framework of a solid-on- 
solid model which is solved numerically. Analysis of the fluc- 
tuation contributions observed for hydrogen adsorption onto 
gold substrates allows to determine the surface tension of the 
free hydrogen film as a function of film thickness. It is found 
to decrease sharply for film thicknesses below seven atomic 
layers. 

Solid surfaces exposed to a gaseous environment at- 
tract the gas molecules by virtue of molecular forces, giv- 
ing rise to a thin film of adsorbed molecules under almost 
any condition. The most important attractive force in- 
volved is the van der Waals force, which is due to the 
quantum mechanical zero point fluctuations of the elec- 
tron shells of the molecules. It is present for all materials 
and extends over distances much larger than a molec- 
ular diameter. It was thus tacitly accepted for a long 
time that the theory of Frenkel, Halsey and Hill (FHH) 
, which is based on the van der Waals force alone, 
accurately describes adsorption isotherms, i.e., the de- 
pendence of the thickness of the adsorbed film as a func- 
tion of the partial pressure of the adsorbed species in the 
surrounding medium. It predicts that the thickness, d, 
of the adsorbed film varies as 
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where a is the van der Waals constant of the system, p 
is the gas pressure, and po is the saturated vapor pres- 
sure of the film material. The power 1/3 on the r.h.s. 
of eq. [I] stems from the fact that according to the Lif- 
shitz theory of dispersion forces, the adsorption energy 
should scale as the inverse third power of the film thick- 
ness 0||. This holds for the film thicknesses to be dis- 
cussed here, for which retardation effects can be safely 
neglected Deviations from FHH theory are to be 

expected only when the temperature is close to the criti- 
cal point of the adsorbing substance, such that the width 



of the free film surface may be comparable to the film 
thickness. 

On the contrary, experimentally determined adsorp- 
tion isotherm revealed considerable deviations from FHH 
theory for virtually all systems studied [[7 -10|. Although 
it has been mentioned before that thermal fluctuations 
may be the cause for the observed discrepancies |l^ , |lT| , 
this possibility has not yet been studied in detail experi- 
mentally. It is the purpose of the present paper to show 
that thermal fluctuations of the free film surface are very 
probably the cause of the above mentioned discrepancies. 
If these fluctuations are taken into account, the adsorp- 
tion isotherms may be completely reconciled with theory. 

As a model system, we have chosen hydrogen films ad- 
sorbed on gold substrates. Films have been prepared in 
a helium flow cryostat at temperatures around the triple 
point, ranging from 10 to 14 K (T 3 = 13.952 K). The 
thickness of the adsorbed films has been determined via 
surface plasmon resonance spectroscopy, as it was done 
successfully before Q . The substrates consisted of ap- 
proximately 50 nm thick gold films prepared by evapora- 
tion onto polished glass. The adsorbed hydrogen shifts 
the surface plasmon resonance of the metal surface, which 
allows to determine the coverage, i.e., amount of mate- 
rial adsorbed per unit area jl2|,[l3| • This can be converted 
into an effective film thickness assuming the film to have 
bulk material density. 

Inspection of the gold substrates by scanning tun- 
nelling microscopy revealed a roughness and a lateral cor- 
relation length similar to what was found before by other 
authors 14 lq]. The surfaces exhibit large crystal facets, 



such that at temperatures below the roughening transi- 
tion of hydrogen (w 9 K flf| ), two dimensional phase 
transitions in adsorbed single layers of hydrogen |L6 17 
could be easily observed. From the substrate topogra- 
phy, it could be estimated by the Kelvin equation |is|] 
that measurements of coverage might lose reliability at 
thicknesses well above 5 nm due to capillary condensa- 
tion effects. Consequently, only films with a thickness up 
to 5 nm were used for analysis. As it will be seen be- 
low, effects of fluctuations are most important at much 
smaller thickness, where roughness effects can be safely 
neglected. 

To start with, we mention that both above and be- 
low the triple point temperature, T3, the thick hydrogen 
film which adsorbs close to saturation is in the liquid 
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phase jlMjill. This is seen qualitatively from the fact 
that if one takes for po(T) in eq. [j] the saturated vapor 
pressure of the liquid phase, which must below T3 be ob- 
tained by extrapolation of the liquid/vapor coexistence 
line, all isotherms fall onto a single master curve. This 
is demonstrated in fig. [I], which shows the data of ten 
isotherms obtained for temperatures ranging from 10.190 
K to 14.067 K, plotted as a function of reduced pressure, 
p/po(T). Since the saturated vapor pressure of the solid 
phase is always below that of the liquid phase, the di- 
vergence of the isotherm is not reached for temperatures 
below T3, such that the film thickness at solid/vapor co- 
existence is finite. This thickness at saturation increases 
with increasing temperature, and diverges for T — > T3. 
This is triple point wetting |Q . 

As one can see in fig. [j], the mutual agreement of the 
isotherms is very good for sufficiently thick films. The 
van der Waals constant is found to be 7±3 K(nm) 3 , which 
is consistent with the theoretical estimate of 8 K(nm) 3 
. However, closer inspection of the data reveals sys- 
tematic discrepancies which are particularly obvious at 
thicknesses below 2 nm, as it can be seen from the inset. 
This is revealed more clearly in fig. ^, where the data of a 
single isotherm (T = 13.964 K) are plotted such that the 
behaviour described by eq. |l| would yield a straight line. 
The characteristic negative curvature of the isotherm in 
this plot demonstrates the deviation, and is similar to 
what has been described before by other authors [||,[l0| . 

The data are very well represented by the solid curve, 
which is obtained from a refined model taking into ac- 
count thermal fluctuations in film thickness p!l| . The 
statistical approach is based on a SOS-model which has 
extensively been applied to multilayer-stepped adsorp- 
tion isotherms ]2l|-p3]|. We assume the substrate to con- 
sist of a square lattice of N adsorption sites i = 1, . . . , N 
and thickness di = aorii of the adsorbed film at the site 
i allowing n, > to be any non-negative integer where 
do denotes the monolayer thickness. The restriction of 
gas molecules to certain lattice sites is well proven for a 
monolayer and seems to be justified for films only few 
layers thick, i.e. the regime which we are focussing on. 
The vapor is considered to be a homogeneous reservoir 
of molecules with chemical potential /i = kBT\og(po/p) 
and the adsorbed molecules are assumed to pile up at 
each site in columns, without forming overhangs or va- 
por bubbles, which is reasonable for thin films and tem- 
peratures well below the critical point. The statistics 
of the film thickness is then given by the partition sum 
Z = J2{ n } ex P I - where the sum runs over all 
configurations rii (i=l,...,N) of the film. The Hamilto- 
nian reads jn| 

N / m \ 

i V P v=l v } <ij> 

(2) 



where z v is the distance of the f-th layer to the substrate. 
The surface tension 7 of the film-vapor interface takes 
into account the molecular interactions within the film, 
where sum runs over nearest-neighbor sites < ij > only. 

Assuming that fluctuations in the film thickness are 
not relevant, one may minimize the energy (Q) yield- 
ing the most probable thickness d given by the FHH- 
isotherm, Eq. ([!]). But the film- vapor interface is always 
undulated due to thermal fluctuations which become im- 
portant for thin films where fluctuations are hindered by 
the substrate. Thus, in order to perform the partition 
sum we apply a mean-field approximation, replacing dj 
in Eq. (0) with its mean value d. One obtains a self- 
consistent equation for the mean thickness, which can be 
solved by standard numerical procedures and can be ap- 
plied to the experimental data in a straightforward man- 
ner (see Ffig. ^) . The solution is the adsorption isotherm 
d([i; a, ao, 7) depending on the Hamaker constant a, the 
monolayer thickness do, and the surface tension 7. Of 
course, 7(d) depends itself on the mean film thickness 
and may be determined experimentally when vapor pres- 
sure p and thickness d is measured. 

This SOS-model reproduces the results of Brunnauer, 
Emmett, and Teller as well as the FHH model in the 
monolayer and thick-film regimes, respectively, but for 
intermediate coverages, a qualitatively different behav- 
ior is found which is governed by thermal fluctuations 
of the film thickness, and is therefore determined by the 
surface tension, 7. As one can see in fig. ^J, the substrate- 
induced hindrance of such fluctuations increases signifi- 
cantly the mean thickness d as compared to the most 
probable thickness d given by the FHH-isotherm. Most 
importantly, the isotherms obtained experimentally are 
well described taking fluctuations into account, without 
any further assumptions. Qualitatively, the effect of the 
fluctuations, i.e., the deviation from the FHH isotherm, 
increases with increasing temperature, as expected. For 
illustration, we have included in fig. |2| what is obtained 
from our model if fluctuations are absent (dashed line). 
This was done by setting 7 — > 00 in the calculations. The 
impact of the fluctuations becomes particularly obvious 
this way. 

Let us now go a step further and use the fluctua- 
tion contributions identified above to determine the sur- 
face tension, 7, of the hydrogen/ vapor interface. This 
turns out to be mathematically difficult on the basis of 
eq. ^. Therefore, we used a somewhat simpler approach 
based on the fluctuation-dissipation theorem, which di- 
rectly yields a more convenient expression for the excess 
chemical potential due to capillary waves j25j, and also 
contains 7 as a parameter. It is less accurate than eq. ^ 
because it neglects the geometric limitation of the fluc- 
tuation amplitude by the finitencss of the film thickness, 
but it suffices for a rough evaluation of 7. 

The result is shown in fig. |3|. While one would naively 
expect 7 to be constant, we see that it is dramatically re- 
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duced for small thickness. The jump at 0.7 nm may well 
indicate the boundary between the solidified part of the 
film close to the wall jl2],[l9) and the liquid. At higher film 
thickness, a strong increase of the surface tension is seen, 
which approaches the bulk value above 25 Angstroms, 
corresponding to about 6.5 molecular layers (one layer is 
3.8 Angstroms thick). We tend to be cautious as to the 
actual amount by which the surface tension of the thin 
films is reduced, because the method used for its deter- 
mination is not as accurate as eq. ^, as mentioned above. 
However, there is no doubt according to our qualitative 
result that the surface tension is considerably less for the 
thin film than for the bulk. It is desirable to develop an 
analysis based on eq. || to check the result displayed in 
fig. I 

In conclusion, we have demonstrated, with hydrogen 
adsorbed on gold as a model system, that adsorption 
isotherms may be completely reconciled with Lifshitz the- 
ory if thermal fluctuations of the free film surface are 
taken into account. This seems to resolve a long standing 
discussion in the field of wetting forces, and also provides 
a means of determining the surface tension of molecularly 
thin liquid films. For hydrogen, we found a substantial 
reduction of the surface tension with respect to the bulk 
value. The exact physical nature of this effect must be 
left to further studies. 
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FIG. 1. Assuming the adsorbed film to be liquid, it is nat- 
ural to use the saturated vapor pressure of the liquid phase in 
the FHH equation. In fact, this makes all isotherms fall onto 
a single master curve, independent of temperature. The slight 
deviations in the flat parts of the curve (around p — 0.5po) 
are significant and may be well explained by thermal fluctu- 
ations. They are shown clearly in the inset, which displays 
isotherms for 10.190 K (squares), 13.939 K (triangles), and 
13.964 K (circles). 
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FIG. 2. A typical isotherm, obtained at T = 13.964 K, 
plotted such that an FHH isotherm according to eq. [I] would 
yield a straight line. The negative curvature is clearly visible. 
The solid curve is a fit of our model. For the dashed curve, the 
surface tension has been set to infinity in the model, in order 
to suppress fluctuation effects. The difference between the 
data and the dashed curve thus demonstrates the significant 
impact of fluctuations. 

FIG. 3. The surface tension of hydrogen as a function of 
film thickness, as obtained from the analysis of the fluctua- 
tion contribution to the equilibrium film thickness. At small 
thickness, the surface tension is dramatically reduced, and 
approaches the bulk value only above 2.5 nm (about 6.5 ML). 
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This figure "FHHplot.jpg" is available in "jpg" format from: 
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This figure "MasterIsotherme.jpg" is available in "jpg" format from: 
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This figure "SurfaceTension.jpg" is available in "jpg" format from: 
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